Statistical Isotropy of CMB Polarization Maps 
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We formulate statistical isotropy of CMB anisotropy maps in its most general form. We also 
present a fast and orientation independent statistical method to determine deviations from sta- 
tistical isotropy in CMB polarization maps. Importance of having statistical tests of departures 
from SI for CMB polarization maps lies not only in interesting theoretical motivations but also in 
testing cleaned CMB polarization maps for observational artifacts such as residuals from polarized 
foreground emission. We propose a generalization of the Bipolar Power Spectrum (BiPS) to polar- 
ization maps. Application to the observed CMB polarization maps will be soon possible after the 
release of WMAP three year data. As a demonstration we show that for E-polarization this test 
can detect breakdown of statistical isotropy due to polarized synchrotron foreground. 



I. INTRODUCTION 

In a very near future we are going to have the first 
"full" sky CMB polarization maps. The wealth of infor- 
mation in the CMB polarization field will enable us to de- 
termine the cosmological parameters and test and char- 
acterize the initial perturbations and infiationary mech- 
anisms with great precision. Cosmological polarized mi- 
crowave radiation in a simply connected universe is ex- 
pected to be statistically isotropic. This is a very im- 
portant feature which allows us to fully describe the field 
by its power spectrum that can have profound theoret- 
ical implications for cosmology. Violation of statistical 
isotropy (SI) in CMB polarization maps is going to be 
very important soon. It can now be tested with CMB 
polarization maps over large sky fraction. Importance of 
having statistical tests of departures from SI for CMB 
polarization maps lies not only in interesting theoretical 
motivations but also in testing the cleaned CMB polariza- 
tion maps for residuals from polarized foreground emis- 
sion. Unlike the foregrounds in temperature anisotropies, 
polarized foreground emissions on large scales. In these 
scales we expect to see the primordial i?-mode due to in- 
flationary gravitational waves at all frequencies. A robust 
discriminator between the primordial polarized radiation 
and polarized foreground emissions is the test of SI. In 
this paper we study statistical isotropy in its most gen- 
eral form based on the Bipolar Power Spectrum (BiPS) 
that was proposed as a measure of SI violation in CMB 
temperature iIi i3, i3|. The BiPS has been apphed to 
check for the SI of CMB temperature maps based on the 
WMAP first year data 0, U- We present a simple for- 
malism that works for all three scalar fields that describe 
CMB temperature and polarization, T, E and B. Then 
we use BiPS as a diagnostic tool to check for departures 
from SI in E and B polarization modes as well as the 
cross terms such as TE. We present an example of ap- 
plying the method on simulated CMB polarization maps 
that include polarized foreground from the synchrotron 
emission in our galaxy. 

The rest of this paper is organized as follows: section II 



is a very brief introduction to polarization and tempera- 
ture anisotropy of CMB and shows how CMB anisotropy 
can be fully described by three scalar fields, T, E and 
B. Section III is dedicated to formulation of statistical 
isotropy in general. Section IV defines an unbiased esti- 
mator for bipolar power spectrum (BiPS) which is shown 
to be a strong tool for testing departures from statistical 
isotropy in a given map. And finally section V describes 
an example of how this method works for a E-polarization 
where statistical isotropy is violated due to large galactic 
foreground from synchrotron emission. We provide some 
useful mathematical relations in the appendix. 



II. CMB ANISOTROPY AND POLARIZATION 
MAPS 

CMB anisotropy is completely described by its tem- 
perature anisotropy, and polarization. Temperature 
anisotropy is a scalar random field, Ar(r;,) = T{fi) — Tq, 
on a 2-dimensional surface of a sphere (the sky), where 
n = (6, (j)) is a unit vector on the sphere and Tq = 
J ^%^T{n) represents the mean temperature of the CMB. 
It is convenient to expand the temperature anisotropy 
field into spherical harmonics, the orthonormal basis on 
the sphere, as 



AT{n) 



l,m 



aimXimin) 



where the complex quantities, aim are given by 



airn = / dririY;*„(r7,)AT(ri). 
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CMB polarization field is described by the Stokes pa- 
rameters, Q{n) and U{n), which depend on the choice a 
local Cartesian patch the coordinate on the sky. One can 
combine these Stokes parameters into two complex quan- 
tities, Q{n) — i U{n) and Q{'n) -\- i U{h) which transform 
like spin-2 fields under rotations of the coordinates by an 
angle tp, 

{Q{n) ± I U{n))' = eT2^^(Q(r^) ± i U{n)). (3) 
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One may thus expand each of them in terms of spin- 
weighted spherical harmonics, ±2^imj 



Q{n)-iU{n) = ^ a2,i„i 2i^m("-) (4) 
Q{n) + iU{n) = Q-2,h« -2Yira{n) 



Im 



Applying spin-lowering (spin-raising) operators 3(3) 
twice on ±2P{n) — Q{n) T * U{n) one can construct two 
spin-zero fields, 



2P{n) 



9 1 -2P{n) 
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For fuUsky maps, the above spin-2 fields can be linearly 
combined to construct two scalar fields 0, 

Ein) = i[9|2P(n)+g|-2P(n)] (6) 

Bin) = i[g^^,p(^)_g|„,p(^,)]. 

Now, expanding these in terms of spherical harmonics, 

E{h)=Y.''frnYim{ny, i?(n)-^af„r,,„(n) (7) 



we get, 



(a2,im + a-2,im) ; 0,fr> 



2i 



{o.2,lm — a-2,im) 



(8) 

Therefore one can characterize CMB anisotropy in the 
sky maps by three scalar random fields: T(fi), E(fi), and 
B{fi) with no loss of information. For cut-sky, E{fi) and 
B{n) mode decomposition is not unique |3jl3- since 
mixing is linear there always exist two linearly indepen- 
dent modes. It is possible to formulate the SI of these 
linear independent modes. Statistical properties of each 
of these fields can be characterized by iV-point correlation 
functions, {X{ni)X{n2) ■ ■ ■ X{nn)). Here the bracket de- 
notes the ensemble average, i.e. an average over all pos- 
sible configurations of the field, and X{n) can be any of 
the r(n), E{h), or B{n) fields. CMB anisotropy is be- 
lieved to be Gaussian [3, ■ Hence the connected part 
of A^-point functions disappears for N > 2. Non-zero 
(even-A^)-point correlation functions can be expressed in 
terms of the 2-point correlation function. As a result, 
a Gaussian distribution is completely described by two- 
point correlation fimctions of A(n), 



C^^'{n,n') = {X{h)X'{h')). 



(9) 



Equivalently, as it is seen from linear relations in eqns. (O 
and lO, for a Gaussian CMB anisotropy, are Gaus- 
sian random variables too. Therefore, the covariance ma- 
trix, {ahnCLp'm') ' fully describes the whole field. Through- 
out this paper we assume Gaussianity to be valid. 



III. STATISTICAL ISOTROPY 

Two point correlations of CMB anisotropy, 
{fii, 712), are two point functions on x 5^, 
and hence can be expanded as 



C'^'^'inu n2) = 



E 

h,i2,e,M 



^afliih^eM^i'^i' "2) 



(10) 



Here ^^^^^^^ are coefficients of the expansion (hereafter 

BipoSH coefficients) and Y^^'j^{ni, 712) are bipolar spher- 
ical harmonics defined by eqn. IjAljl . Bipolar spherical 
harmonics form an orthonormal basis on x and 
transform in the same manner as the spherical harmonic 
function with £, M with respect to rotations ^3- We can 
inverse-transform C-^^ {fii, 712) in eqn. IjlUI) to get the 
coefficients of expansion, ^^fj;^;^' multiplying both 

sides of ean. lflOl) by Y^i^p {hi, 712) and integrating over 
all angles. Then the orthonormality of bipolar harmon- 
ics, eqn. (|A2|) . implies that 



A 



XX' 

£M\hl2 



dfin.j da^.c^"" (ni, h2)Y;i'/^hi,h2). 

(11) 

The above expression and the fact that C'^^{hi, 712) is 
symmetric under the exchange of fii and n2 lead to the 
following symmetries of ^^^j^j^ 

aXX _ (_.Uh+h-L) aXX (,r,^ 

aXX _ t.xx X I. — ^ -i 

^e.M\U — ^tM\ll 0^,2fe-l, K — i, /, O, • • • . 

It has been shown |n| that Bipolar Spherical Harmonic 
(BipoSH) coefficients, yl^^j^^^, are in fact linear combi- 
nations of off-diagonal elements of the covariance matrix. 



A 



eMlhh ~ E ( ^)™^^limil2-m2- (1^) 



where Cf*^ , are Clebsch-Gordan coefficients. This 
clearly shows that ^^^^j^;, completely represent the in- 
formation of the covariance matrix. When statistical 
isotropy holds, it is guaranteed that the covariance ma- 
trix is diagonal, 

(«L«rw) = Cf'^'Si.S^^, (14) 

and hence the angular power spectra carry all informa- 
tion of the field. Substituting this into eqn. Ill3() gives 



^fmw = i-^yc^'' (2^ + 1)'/' Sw S,o Sm,. (15) 

The above expression tells us that when statistical 

isotropy holds, all BipoSH coefficients, ^^^Vfjiz'' 

except those with £ = 0,M = which are equal to the 
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angular power spectra up to a (— 1)'(2Z + 1)^/^ factor. Bi- 
poSH expansion is the most general way of studying two 
point correlation functions of CMB anisotropy. The well 
known angular power spectrum, Ci is in fact a subset of 
the corresponding BipoSH coefficients, 



XX' 



(-1)' 



-.A, 



XX' 



(16) 



Therefore to test a CMB map for statistical isotropy, it is 
enough to compute the BipoSH coefficients for the maps 
and check for nonzero BipoSH coefhcients. Every statis- 
tically significant deviation of BipoSH coefficients from 
zero would mean violation of statistical isotropy. In the 
next section we discuss this in more details. 



IV. UNBIASED ESTIMATOR 

In statistics, an estimator is a function of the known 
data that is used to estimate an observable quantity. An 
estimate is the result of the actual application of the func- 
tion to a particular set of data. Different estimators may 
be defined for a given observable. The above theory can 
be used to construct an estimator for measuring BipoSH 
coefficients from a given CMB map as. 



aXX' _ 
^iM\a' — / . 



(17) 



Non-zero components of BiPS imply break down of sta- 
tistical isotropy, and this introduces BiPS as a measure 
of statistical isotropy, 

STATISTICAL ISOTROPY =^ ni = ye ^ 0. 

(21) 

It is important to note that although BiPS is quartic in 
aim, it is designed to detect SI violation and not non- 
Gaussianity [1 H H [H [H . 
BiPS is given by 



An un-biased estimator of 



r^xx 



El A^^' 
\^eM\ii' 



ll'M 



53 



XX' 



(22) 



where 23^'^ is the bias related to the SI part of the map 
and given by the angular power spectrum, Ci, 
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XX' 



{kf) 



(2£+l)^ ^ Wi.Wi, 

h /2 = K-ii 



riXXr<x'x' 



(•23) 



The above expression for 25^ is obtained by assum- 
ing Gaussian statistics of the temperature fluctuations 
[ll lllj|. Note, the estimator kt is unbiased, only for SI 
correlation. In that case, ensemble average of ki is same 
as its true value which is zero for i 0, i.e., (ki) = 0. 



where Wi is the Legendre transform of the window an 
isotropic smoothing function that can be applied to the 
data. The ensemble average of this estimator is given by, 



V. EXAMPLE: POLARIZED SYNCHROTRON 
CONTAMINATION 
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(18) 



which is its true value. Akin to the well known quadratic 
estimator Ci = ^jqrf J2m Wim]"^ C*;, the above estima- 
tor is an unbiased estimator of BipoSH coefficient. How- 
ever it is impossible to measure all Afj^^^, individually 
because of cosmic variance. Combining BipoSH coeffi- 
cients helps to reduce the cosmic variance. Among the 
several possible combinations of BipoSH coefficients, the 
Bipolar Power Spectrum (BiPS) has proved to be a useful 
tool with interesting features. BiPS of CMB anisotropy 
is defined as a quadratic contraction of the BipoSH coef- 
ficients 



i^xx' 



E 

1,1', M 



A 



XX' 

iM\ll' 



> 0. 



(19) 



Some interesting properties of BiPS are as follows: it is 
orientation independent, i.e. invariant under rotations of 
the sky. For models in which statistical isotropy is valid, 
BipoSH coefficients are given by eqn 116|l . and therefore 
lead to a null BiPS, i.e. = for every l> Q, 



As an example of how one can detect deviations from 
statistical isotropy in CMB polarization maps, we make 
statistically anisotropic polarization maps and estimate 
the BiPS from them. This can be done in many different 
ways but here we choose a simple method which results in 
severe violation of SI and therefore is good for a demon- 
stration of the method. We add polarized synchrotron 
emission template to the background CMB polarization 
map. Polarized synchrotron template (30 GHz) is made 
using Planck Simulator which uses the model by |Q| , 
i.e. the polarization degree is a function of the intensity 
spectral index while polarization angles are derived from 
a gaussian distribution. Here we restrict our attention to 
E mode polarization only. It is obvious that everything 
can be done in the same way for B mode as well. The 
estimator will then be 



aEE 



E E 



c. 



IM 



(24) 



where af^ are the spherical harmonic transform of the 
background CMB polarization map plus the polarized 
synchrotron radiation, 



^xx' 



(20) 



+ a 



Esy 

Im 



(25) 
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FIG. 1: CMB polarization maps with no foregrounds are 
statistically isotropic and have null bipolar power spectrum 
(top). Adding polarized synchrotron emission violates sta- 
tistical isotropy and results in a detectable non-zero bipolar 
power spectrum (bottom). Red dots show the BiPS after bias 
subtraction and Blue lines show the I — cr of the cosmic vari- 
ance. 



FIG. 2: Adding 5% of the polarized synchrotron emission just 
violates statistical isotropy (top). Adding 7.5% of the polar- 
ized synchrotron emission clearly violates statistical isotropy 
and results in a detectable non-zero bipolar power spectrum 
(bottom). Red dots show the BiPS after bias subtraction and 
Blue lines show the 1 — ct of the cosmic variance. 



and Wi is a isotropic filter that allows us to target angular 
scales of interest by filtering out power on other scales. 

We simulate 1000 statistically isotropic CMB polariza- 
tion maps, add the synchrotron template to each of them 
and compute the BiPS for them using the estimators of 
eqns. H24|l and (|22|l . Filters that we use here can be di- 
vided into two categories: low pass Gaussian filters 

<=iV^exp|-(|±i)'} (26) 

that cuts power on scales, I < Is and band pass filters of 
the form 

(27) 

that retains power on scales Is < I < h, where Jo is 
the spherical bessel function and N'~^ and are nor- 
malization constants chosen such that, ^2i(;+iJ'' ^ 
i.e., unit rms for unit flat band angular power spectrum, 
l{l + 1) C^^' = 27r. 

Results of this computation are shown in Fig. ^ and 
Fig. 121 We see that CMB polarization maps with no fore- 
grounds are statistically isotropic and have null bipolar 
power spectrum. Adding polarized synchrotron emission 
violates statistical isotropy at large angular scales and 
results in a detectable non-zero BiPS. Retaining only 5% 



of the polarized synchrotron emission just violates statis- 
tical isotropy at the threshold of 1 — cr. At 7.5% of the po- 
larized synchrotron emission clearly shows the violation 
of statistical isotropy and results in a sharply detectable 
non-zero bipolar power spectrum at K4. 

We should emphasize that this is simply an example 
to demonstrate how violation of statistical isotropy can 
be quantified in CMB polarization maps. In reality, we 
usually expect to deal with cleaned polarized maps which 
would contain some residuals that have different angular 
structure. The signal would be much weaker and also 
have different BiPS characteristics. Hunting tiny residu- 
als from foregrounds in maps of temperature anisotropy 
using statistical isotropy has been studied 15] and sim- 
ilar strategy can be applied to polarization maps when 
they are available. In addition, other observational ar- 
tifacts such as anisotropic noise or incomplete (masked) 
sky can also cause violation of statistical isotropy in a 
polarization map. In the latter case, the incomplete sky 
coverage immediately induces a contamination of E-mode 
of polarization by its B-mode and vice-versa. Then the 
modified temperature and polarization fields is related to 
their actual values of full sky coverage by a window ma- 
trix 1^ 01 whose elements are basically window functions 
for temperature and polarization in harmonic space. It 
can be shown that the estimated BipoSH coefficients are 
in fact linear combinations of that for fuUsky CMB maps. 
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Here bold-faced A^j^^i^^, and A^j^^^^^^^ are the col- 
umn matrices corresponding to estimated and true Bi- 
poSH coefficients respectively, for the auto and cross- 
correlations (TT, TE, TB, ET, EE, EB, BT, BE, BB) of 
temperature anisotropy and polarization. The elements 
of the matrix N^^"j^;^ depend on Clebsch- Gordan coef- 
ficients and window functions in harmonic space. Hence, 
the true BipoSH coefficients can be estimated from the 
pseudo-BipoSH coefficients by inverting the above equa- 
tion. We defer this to a future publication, a SI analy- 
sis of CMB polarization when this effect is important. 
However, we have verified using simulations that the 
BiPS of polarization maps is insensitive to breakdown 
of SI due to galactic cut when it is filtered at low-1 us- 
ing, W^{h = 10, /* = 0) and Wj^ih = 20, k = 10) in 
eans. lf2B|) & (j?7jl . (This is consistent with the result for 
cut-sky CMB temperature maps discussed in the paper 
[Tl|.') As a result, the BiPS signature of the polarized 
galctic foregrounds presented here would not change if 
the maps are masked by a galactic cut. CMB polarization 
maps filtered with windows peaked at higher multipoles 
(eg., Wf{ls = 90, k = 80)) do reflect the SI violation aris- 
ing from a galactic cut. The complications of quantifying 
statistical isotropy in cut-sky polarization CMB maps is 
formally encoded by eqn. I|28(l but its implementation is a 
challenging task which is currently under progress. (The 
effects can also be estimated through extensive simula- 
tions.) 

VI. SUMMARY 

We present a novel approach to quantify the violation 
of statistical isotropy in CMB polarization maps for the 
first time. We present a fast and orientation independent 
method which allows for a general test of isotropy using 
Bipolar Power Spectrum. This method has been previ- 
ously applied to the temperature anisotropy maps and 



many various aspects of that are well studied in details. 
In this paper we extend BiPS to the CMB polarization 
maps and by present a working example to demonstrate 
its potential. 
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APPENDIX A: USEFUL MATHEMATICAL 
RELATIONS 

Bipolar spherical harmonics form an orthonormal basis 
of X S'^ and are defined as 

^W'('^l>2) = '^fimii2m2^'i™i (^1)^2^2(^2), 

mim2 

(Al) 

in which Cf^^i^^^ are Clebsch- Gordan coefficients. 
Clebsch-Gordan coefficients are non-zero only if triangu- 
larity relation holds, {hh^}, and M = mi -|- m2- Where 
the 3j-symbol {abc} is defined by 

r Li_/l if a + 6 + c is integer and |a — 6| < c < {a + b), 
^ ' [0 otherwise, 

Orthonormality of bipolar spherical harmonics 

(A2) 
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